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PemelicsoChiFITON CF TEE FROBIES 


ier mOmmner Stas eh VGRhCUyiamic E€CUadticns Gdescribrin, 
Saeemechic crlow cver a cone in ait has cFeen well kKrcwo Sirce 
moe. 1530's. Unt recently, the problem of déscrikinge tae 
flow cveér a cone in sater has been limited to scluticrs of 
the sukscnric case. PElMeGai Caley lations were ligited 
subscric ilcw crcecausée researchers relieved that surersoric 
flow in sater was not reasiltle for normai vessels (Such as 2a 
shap). PormmOmecliakhy VEssSeuseiMewatOL, 1t 1S cCertairiv true 
Meee sULeELSscnic flow fast that vessel is hijgkiy imprerabie. 
HOWEVEL, fio MeL LOR eOhwtne = Nera get Croft am expicsive 
shaped-clrarce rired trderwater is supersonic. 

This thesis presents a Cecmputéer program which calculates 
M@emenyaorLcdyzramic fics Fast a cone im either water cr ai 
tee cr Sstrersonic cCenditions. Tne program utilizes the 
methcds developed by frevicus researchers £cr calculatin 
the superscnic flow in air and whicn have been suitarl 
MOdified tc describe tne ccnditions in the water. ree | 
[Pemmededtichs 1nclude utilizing the moditied Tait eguaticaz, 
Meee i> £06 “tnermnal" cr “thermodynamic” e, uation cf state 
meme water, to describe the  ~nysical stare of tne sater 
Batner trar tae ferfect gas law used for air. 

Poem lchicene tLe Gone GiMer in the jet from an expio- 
Pave Shaped-charge is blunt-nosed rather than ccrical. 
HOWeEVEL, Serronmmer se the ccenreal case iS a prelifinary 
requirement to soluticn of the actual blunt-nosed fretlen. 
The sciuticn to the conical flow case will serve as an 
meee cot test prLOgrEan ror the solution to the plurt-rcsed 


Frobler. This thesis presents a solution to the ccrical 
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Tke ccomputer precram presented if this werk was "Cie 
nally developed in the BASIC computer language usirg a 
Hewlett-Fackard HP-67 cCenouer That program is tne Lasis 
for this thesis. It was desired to translate the frecgran 
into a kigter-order ccmputer language Lok (~€xecu tics vere 
larse, wainirane ccoputer systea. Innis translaticr Was 
desired in crder to gake the precogram more readily acc¢ssitie 
to a wider pody of researchers and if order to spfeec the 
executicr time of the prcgrag. In tnis tnesis, the 


Iclicwitsg gcals have teen accomplished: 


« Successfully trarslate the fErogram tren SAStC toe 
higker~crder lancuage. This jgoal was net by utiliaie 
FCFIFAN as the tign-level languaye of choice. anile 
FCRTRAN has gany drawkhacks as a high-order language, it 
is still widely used in the scientific ccinunmae 
FCRISAN was USed, therezore, SCmten ee ieee emer 

e 


useful to as wide an aréa cf researchers as => 


e Felicw nodern picgramming ~practices in the design wae 
jauplerertaticn of tne pregram. AS Beltore, tre eheiwec ves 
FCRIFAN as tke high-levei language maxes this gcal some- 
wodt UCLeMesieer Tonite However, Many computer scientists 
have demonstrated that structured froyramming fractices 
can fe achieved tSing FORTRAN. To the largest extent 
gEcssatle structtred programminy fractices have ce€een 


utilazed. 


e FIesent a "user-friendly", well-documented prograr. In 


this regard, literal use Of ~eomments ceocus, in Vee 


ar = -«- - t= oe i ee NN NS oe a Se ee oe = Seo ate 
Paw Go S toe ee, eee ee el (aLleaatvae eed = od ace WoC ug fae! 
epee Oa eno tom—e at Gn SES Olstrate tae Losic ci ts 
Weeccler <2 lewimee eeca- i eeeciceccene OT, Wewro Lidgite 


ORCC ORM nee ene USer, the user's restcnses ar 


verified befcre tre prograd executes. 
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Il. FUNDIN ENTIAL soue ame 


A. FACKCGCECUAND 


Tke first step necessary toe descrise surerscric flow 
cver a cecne igs to calculate tne thermodynamic freperties 
across a Shceck wave. The basic research into the Chancemam 
thergmcdyranic preperties across a shock wave in water was 
extersively concucted and reported upon in Uncerwavem 
Explesicr Researca [hef. 1] during and jWsSt after wer) Gee 
Hi Te Toe b-est Summarization cf tnese works can ce fcurd in 
Cole [kef. 2]. 

The primary source used as reference for the calculation 
ci tke Eydredynamic froperties of Séa water at Ee Li creem 
asacck wave is the scerk of Richardson, Aronus, and Halverson 
(Reis. age they utilized graphical technigues, ~ whieaweeee 
rather crude ang tecicus, te calcuiate the tnerredytraire 
data needed te describe the conditions of the sea water. 
Fuhs [Ref. 4] used tle work ci Ricnardson, ¢€t al., (hoi eee 
to develoz; a ccMputer program for the HAPYICY anand neue 
Ccalculatcr shich efficientiy calculated the same therigeys 
banhic Tero eer eave Funs' jERef. 4] programs provided Waa 
Fasis fcr tae FORTRAN subroutines, used ifn this work, wages 
calculate tle same thkermodynatic froperties such as grease 
sure, tetlperature, arc density. 

Ibe funcamental ecuations used to calculate the tnerimo- 
dynamic charges shick cccur across an oblicue Shocx If¢crte 
air héeve Leen well krcewn fcr decades. These equaticns are 
excerticrally well described in Kinney and Granam [kef. 5] 
and ifcrm tke basis fcr the initiait Calculators icra 


superscnic flow fast a cone in aic. 
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Pomme ac lait secMaticn is the €eiation mucst orten usei to 
Bescrice the tsaermcdyranic State of the water. The Iecisied 
Tait e€guaticn can be used te describe €ither pure water or 
sea water. The forz cf toe rodified Tait eguaticn used in 


Mmeomchesis WaS LaKkel treme acit [kef. 63], whe nas Ccotintedl 
to perfcrm research ir underwater explosion fhenozena. AL 
excellent discussicn c£ the mcdified Tait eguation ard tow 
at Can ke utilized is contained in Rowiinson [Ref. 7°. 

After tre théermcecynamic ;rcperties cl the water (cr air 
as tke case nay be) cn the dcwnstream side of the shccx 
Meme naye Leen calcusated,, an iteration metnod, utilizing 
Mmm——eritciatea Grapaical CcnsStruction first eveicpec cy 


meeemain, | kei. 8}, is used tc pregress from the shcck frent 


momeene cche surface. Tne equations needed for use fy this 
mecraticr Wethod are fuliy déscriced, LOD Tale; See cera S 
(Rer. 6&6}, who descrikes their development and use. The 


Metoocs which aprly tc tne itlcow past a cone if air can, witn 
meme ceCSSary Changes made for the dirrerences in tre tner- 
modynarics of the tso fluids, be used to calculate the 
Semmecl SuUcreErSOniIC <flow in water. These calculaticns icra 
Metco oO the. Fain pact cl tne rCrRTrAN pregraz which 
<Ol ea os iPoneGteeatesmOtLedm nNele that Snapiro alsc fecirts 
Mier Ge ;LeEReCering werk Gi Tayior and Haccoll in tre 193C's 
and i$40's on the méthods ci sciution of the flcw over a 


meme If alr problem [Fef. 8}. 


eee cENEEAL EQUATIONS 


in tte develcpment of the €guations which follcw, it is 
assumec that these eéquations can be validly used te calcu- 
late the flow over tke cone in either air Or water. The 
Seiieerenhe were ~~ —Erimagily developed by Shapiro (Ref. 2° for 
mugw 21 ait. 


JfFLEI Ca@_couazti Le the CeLEI ICL bells 26a eS eee 
as a result of tre Fhccx crent, it is 4SSuULeQ hace eee Glee 
rrorerti¢es aiil remaiz CorSstant On Ataganpar, cone= £2 Vie 
ccmaocr vertex. By this aSsumrtion, the tlow past (le ieee 
can fe calculated. Tne ilcw geometry is illustrétec by 
Ligdubes2. | sand wiagur cee. which snow a typical streazline 


and its izace in the Ecdograch [lane. 

As discussed therceugniy by Shapiro { ker. dj, there ee 
disccrtirucus change ihn both direction dnd velocity dcreee 
the cklaicue shock frent. Pcints ? and 2 {see figure z. ema 
figure 2.2) lie, tkeéerefore, Ch a Common SniccxK EClaxr “aamem 
CEIGinates at poms Between point 2 and the Pointe 
wnich is cn the cone surface only at inkinite distaneee 
there 18 a YFregicn of comical ilow where the stream ~rezgeg. 
tieS vary ccntinuously. The velocity vector to the “pemine 
's' arn tke nodograpr plane defines waat is called tie wea 
semi-vertex angle witk the centeriine axis. In the wetacds 
Whauem  f Cauley the ccne semi-vertex angle 1S tne variatle 
which is ultimately ceternmined. Further, as pointed curiam 
Shapire [Ref. 8], since all streamlines in the ficw exreri- 
ence tre sage entropy jump acrcess the shock £Eront, the fie 
tetween tre shock ficnht and the cone surface 1s Ecth isén- 
tho Pic ala 2Er © cee acme 

Tke second-crder diicterential eyuation which éctually 
cescrifes tte frlow cf tne fluid fast the’ cone in air is 
Lusty “€eveicped Dy Stasi pom nei - sone Shapiro a2o0tes that 
there are tno methods commonly used to solve this €Cua tue 
Coe mwethcé, developed first by Taylor and Maccoll | Rete 
perfcrms a rumerical integraticn of the equation. The seéeccnd 
methcd, which intecrates the equation using a  Grarpnueas 
constructicn method, was first developed by Buséemarn. MTunae 
Flrogragd developed in this toesis utilizes Buscmann = tetnedw 
rcdified fcr performarce on a nodern high-speed computer, to 
ferficrm the integraticn of the full second-crder differen- 


tial e€gquaticn. 
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Figure 2. 1 Shcck Cone and Typical Streamline. 
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Mee SaGel ceoraltat=es £ aba . séve Teen 313525 «itr tne 
memeees Cla ifG velocity com;crents ies lou - foec, fel ot ls 
eee) Go eetecmmvetiadr, <GlatlOns which are used in tne 
computer ~rcogram are fresented in tais tresis. Pee bolled 
@ecivaticn cr the equaticns can ge weeer naw es. Saaciro 
(Ref. €1], and neti net ce cfepeated nere. The aomenciature 
Geeao aim the equations develcped in tals cnapter is detéeiled 


me darlie I. iene cl i mmew ebm mOocCern ChOUSRE, toe S=K=S 
{meter—-kilcoram-secordg) unit system has been used threusgnacut 


Mepeeetresl=S except icr occasicnai lapses during the sat 


ap) 
ry 


calculaticrs when pressures are referred tc in urits of 
kilcfrars. 


meee ncmaeORCtTEY OL Preure 2.3, 1t Can cé SEEN teat: 


rs 

a 

ond 
ee? 


ue = Vccs(w- 9) and Ve Sn yt ) ( 


Meetene cévelopment ci tne actual second-crdéer dirtierertia. 
Pepemeren, =sfapite f Refi. 8] Enews that, due to the ccndition 
Smee roOtaticnality, tke follcwing relation must ce true: 
dV 
bs Ga 


Meeeeectetiatahg Hon. z.1] with ressect to w , cne cktairs: 


ic d9 . a dV (es } 
oe ci = a au “gon 8 
a vac Au Sih 8) + aR Coe ) 
and 
dV 
wo dg ; ie a er oe 
a * -V(l - =) cos(w - 6) ae aos) ( ) 
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Shock Cone 





@M 
eo 
fas! 
ee 
a 
> 
— 
w 
= 
@ 
O 
os 
oe 
= 
~ 
See 
ee 
a 





a ee ee ee eee eee ee ee ee ee ee ee ee eee ee eee ee eee ee ee se ee ee eee eee eee ee es se eee ee ee eS nd 


Cone and Flow Geometry. 
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Piss ltuticn Of this restizt inte f£jn 2.3 yields: 


ah (2.6) 
dé dw 


dw  Vtan(w - 8) 


PeeweorO f[Ret. 8] shows that the e¢gjuation yoverniny the flew 


Mme CCLicdal Legior is: 








dV 
- 4) Z Z 2 
es ty (av + V cotw + ——)(V - Vv -V-) 
Z is (33 da max 4) is 
a) 
dV 
dv 
rr a 
r dw () yV 
dw “ @ 
—- . e . * d8 —— - oe e 
Elimaratings aa teomemegezes and SUbSTIiItuting the exrres- 
: dV 
5) = f c O) an | me nS 7 5 
earoms ICL vo (ere se eweoOo Oy een = 2.1, aged and 


V 
feo iptc EGn 2.7 gives: 


Sono 
dv _ Vsin(w - 8) sinw 
ee ean near) ra) 
8 2 - 
Coe Eee. -V) 


Mermmmnating "“R* as tre radius Of curvature of toe pncdecraph 


osem@edatlire, che cbtairs: 


dV > 
_,__<:eeee IR et See (2.9 
s sin(w - 8)dw 


1& 


OL ¢. <egaid ves tae. 


Irv) 
é 
iO) 
s 
jo 
ci 
2 
tr 
car 
ri 
NR) 
e 
01) 
% 


ti 


s_- 
ces 


CCTLi ice se 


sing 








V Sone (Zs 10) 
oe ye ae 
29 Sin (4 = oo 
Gee pat) 
max 
Eut, tre energy equation asserts tEnac: 
Z k - lL Zz 2 a 
= ~ Le 11 
———_) (an. (Z. 11) 
PERCE Ee Lei e 
sind 
R = V sind 
ve ee ie) (2. 
| a aceesenea 
2 
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where c is the lccal ¢feed cof sound in the fluid. 

Ecuaticn 2.1z 15S the pasis £06 the calculaticn Come 
superscric flow cver the solid cone in ¢€ither a1 CE wae 
The crapfhical integration method invented by Buserann is 
adegtately explained in Shapire [Ref. 8] and need nct be 
repeated kere. Essertially, the computer proyram given in 
this werk automates the Eusemann graphical integratie 


metuce £c¢r calculatire the ccnée semi-Vvertex angie: 
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The eEgtaticn of state ick ain is Specitied 2yeee 
perfect gas law {under the asSsSumeticn, that 2s, that. tae 
rehaves as a perfect gas). This law is quite elegant ang 
Simple ard allows easy Ganizulation to obtain various Guan- 
tl ti ec. The form cf the perfect gas law mcst often used in 


the calculations of tktis thesis is: 


p = of8T (Zeley 
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MeeevetSedi gaS COrstart, joes 
A a 14 
R 7 tis } 
a 
Memecagudticn 2-147, M, is the geclecularc weight or the air. 


Ir air, the cnarge in the thermodadynagzic properties of 
the cas as it cresses the shcck front are easily calculated, 
as shcwnp in Kinney ard Grahag [Ref. 5]. SiGe es MOL t sri uid 


Gynazics teitbooks iliustrate the develo,ment of the €egua- 


~ 


tions inicn follew, it is net necessary to derive €M nere. 


Te 
AS méerticned previously, Kinney and Graham [ Ref. 5] 
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Mmeeptaicntally lucid explanations and derivations. fhe frin- 


fly 


Ciple ecuations used te calculate the thermodynamic ckances 


which c¢ccur across tke shock Iront in air are as follews: 


en ee Bc 


tan(S - oc) _ 
f tan® 


ess 


5 
(k + 1) M, sin 8 


rm PIir tY 


FQWaticn z.15 if used to iteratively determine the aefiec- 
men ancie @. All cther guantities in tnis equaticr are 
Meee {lec 8 1S tke Shock argle and a € < 
twee DURCeET, bcth ci which are input varameters tc t 


mmecgidi ict the air céelculaticns). 
2 2 

) 2+ (k - 1) My Sin 6 

a ee Zeeetles) 


kM rece - (k = 1) 


[M., 


Bavinc determined @G from e€quaticn 2.15, eguaticn 2.16 is 
used tc determine the Mach nurker on the dcwnstream sice or 
suemeetCcck front. 
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Figure 2.4 Gecuetry cf the Obligue Shock Front. 
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tf) 


meaner IGCEe COmplicated. Tne cSt commonly usea state Eecua- 
tion for water is knewn as the modified Tait equation, skich 


mie Lée written as; 


p = B(S) (CS) a (2021) 


waeeretc Es) iS a Slowly varying function of entropy aicne, 
Memeetresamately a ccrstant equai to 7.15, and ¢€ is the 
value ci the density for zero pressure. The above ¢€cguetion 
Mamet ctorei: | Ret. 6), Dut, in different forts, it is also 
mmm@mm@treac GCUIte extensively in Cole [Ref. 2], Rickarcscni, 
fee al., | hef. 3) ard ROWLlinscn [Ref. 7]. AS= acter ered ih 
Mmemerrt reduction to this work, the moditied Tait eguaticn is 
meme therfai” equaticn of state for water. COE itch: ecay 


shows that the modified Tait E€guation is of the forn 


eo (2.22) 
p = Bs) [ESS = ° 
Grpeeit Siiifier terms, 
kf = plv,T) (eae 2) 


peee tie full modified Tait eguation welates the three 
SiemeGedyralic Quantities of £, v, and T, it is calliec the 


"theracdj;nagic™ cr "thermal" equation or state. 
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Sitilariy, if Ccrt<raést tec ait, een Gas eed a eee eee 


‘ 


wnich cccur aS 4 result o£ the SaCccK DECCESS Ih Water Capa 
be easily delineated [fy siacle ecuations as in tne air case. 
BOWEVEL, a Sinelification can ce made in the water case 
kecauseé, uriixe in tre air case, tne pressure jump acrcss 
the stock in water is so very large. In alr, [L©€sstre 
changes acress the shccxk trent ch the crder cf 1 tc 2s cee 
are ccnsidéred large {at least for chemical expiosSicrs) 2 ie 
contrast, as pointed cut by kichardson, et. al., {het lee 
the Pressure JUMP acrcSs a Sheek if water 15 ch ehe rca 
kKilclars €¢ tensor ilobas <4 Thererore, tne calculaticns 
can ke simplified Ey ¢fecifyirg the pressure on the Comme 
Stredi Side “ORMPBewe i cee: cent This 1S vaiid since tue 
upstream [ressure is sc small in compariscn to the urstréan 
dynanic [Eressure 01 vy 72 and in comparison to Ect: ieane 
downstream fressure ard dynagzic pressure. Tne Specificaticn 


cr the dcwnstreaz Pressure 1S accomplished, in tne poegue 


Kh 


Ci -€b1S store sise by aliowigg the user to input a "~ress tee 
multagplication factor (MFACI)", woich converts the ~rece ee 
upstream tc a pressure downstream at point 2 which is Gives 


ays 


BE = ie & MEACT xX 1Coees © (Z. 24) 
Tne factcr 1000.90 14r egquaticn 2.24 Converts the bricht sae 
side cf the egquaticr from a fFressure in Dars to a fresstre 
inykKi Jobers + AS tan eva cy F, is 1 bar and MFACT is 5.0, 
the pressure at foint 2 downstream will be 5 kilobars. 

It can be snown that the Simpiirtying assumeticn made 
above is exrtirely vaiid by ccnsidering the momentum eguaticn 


for steaay frictionless flow alcng a streamline. The 


24 


See -Girelrsicual qoneLtug eyuaticn states: 
12 ay ea) 
_ - o (2 _ D a = ce 2 Deg = 
ET Bt 72 
Mea necte that typical values for 9, and ¥, Pomme as paocc ler 


under ccrsiceration are as fcllicws: 
©, = 1000 kog/a and Vo = 1500 m/s 


Ther¢icre, vi /2 = ste 107 Based te OF auout lt. | Kilo- 


~ 1 
Mieo. i2yrical values for Pros tne upstcean pressure, are cn 
mremenraer Ci | te perhaps a few tens of cars. Thus, it car 
ke safely assumed that F, can be ignored ccmpared tc tne 
tpstreag dynamic pressure. 

PaveanGe sPeciLieaq Luis ~ELressure, tne program cr tris 
thesis utilizes the subroutines developed ry Fuhs [fFef. 4] 
to calculate the density and velocity components at a forint 
Mist CCwrstream cz tke shock frent. Since the gecretry of 
the velccity flew across tke shock fI[ront is the sare in 
€1lther air cr water (see figure 2.4), phe-velocity CcHnr~e— 
rents determined in tke first step can Le used to "Junge pack 
across" tne shock rcrert to tne upstream side where tie free- 
stréeag Mach number 1 can fre determined. 

After determininc the conditions on the downstreaz side 
Meee STOCK EEont, the program iterates aiong tne strean- 
dines tc determine tre cone séni-vertex angle. This itera- 
tion is performed using the Busemann graphical construction 
as for tke air case. The majer difference Létween tke water 
and air cases, after the initial shocked conditicns nave 
Feen deteruined, Pewenat the local speed Ceo eer d sen a 
streaglire gust be dé€termined iteratively, for the satexz 
case, ky uSing Fuhs' [Ref. 4] subroutines tc deterrine the 


density and pressure at each foint. Knowing the pressure 


and cersaity, the Sre€ed Of ECUTA in the water Can ce Cala 
Lac eG Glee 


n 


NO) 
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Z ie es 
6 Gee =) 


a = 


which is derived fren the modified Tait €quation Yertpewa 


defiriticznz cr the speed Gk Semred as: 





eZ Ip | (Ze 27) 
ec = 
Jp ‘Ss 
These equations, Ccukined with the eguations develcreaqmias 


secticn E— cf this chapter, constitute althe equduaew 


needed tc ccmpletely sclve the supersonic flow over a ccne. 
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Nomenclature 


Specd seme scurd 

Bic teEe of Specific neats ao 
(are Gem Uete9) 2 ie 

NeleeGular weight Of aia 
Pressure 

potemwsomaia Sscherical»cecrdinates 


Pw eemonecULVAtite or Cac 
iedodeamm Streamline 


Temperacure 
Velocity 
fare ve LeClmumincl dddabatic flow 


Angle the stréamline makes with the shock | 
Elenemercame Tiommme wanGleo) 28 air Caiculaticns 


Ccne semi-vertex ancle 

Ficw directicn 

Mass density 

Shc Capargic 

AmeGlembieseterical cccrdinates 
Universal gas constant 


eromgelesed CCudLtich WpStream oF 
Scimeae, SiOCK Cl a ShCCK Front 


Sa wees Secondlrign downstream oT 
CGimeeleshcck Cr aesmeck Tropi 


Signifties ccnditions at cone surface 
SlciieerememmeencCOneCHt im redirection 


Si Chimeeese a —cCMrOnent in w-direction 
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Be SERCCEAE OO Gure 


As discussed in Crapter 1, the computer prograrc devel- 
cped fcr this thesis saS written in the FORTEAN Ercclranaee 
languace. Sir teegumed PLCIURARDL DG aed eee have Leen 
folicwed throughout in that the program is dividedc ages 
ElocKs (Cr Tea aeics,) c£ codé, eacna of which 1s designed 
pericrm a single calculation sequence. Many or the gccules 
have teen pflaced inline rather than dDeing wlitter as sépa- 
rate stkrcutines or ftnctions. This was dene DLidarily Gem 
€ase cf tse and understandirg. In addition, however, mest 
modules were placed inline because tne parameter 11st 
excharced Letween the module and tne main frograz xculd nave 
reen excessively ilcac ctherwise. Regardless of whether the 
code is included inline or as a separate module, the ~regran 
waS wratter in a manner wnich ensured that the "side- 
effects" problems discussed Ey MacLennan [Ref. 9] did ieee 
cccur. In addition, the use of FORTRAN's CCHAOW cCcr= foe 
has Leen studiousiy avoided tc prevent tke aliasing frecrien 
Gascusséed Fy Naclenian rere ae A module hierarchy; Chae 
whick enows the maztcr modulés of the fprocram ana tareuee 
INtELCCKIrEections, 1s Clven ds fee 

Ike flcw logic cf the fain computer program is ieee 
tratec Fy figures A.1 through A? of Wopendax A. Ihese Tite 
charts shew that the initial sart of Che Vrogran (uf): oe eee 
60) 1s used tc initialize certain key variables and to 
gather irput values from the frogram user. Note that, am 
stated in the objectives for this projram, all user in suc 
verified tc ensure that the values entered fall within 


Frescriked limits. In crder to ensure that incut values are 
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Tpecsge ic ails, the SuUrrOuUtiCe CGAasiNe= 2s Devore ce Je. 
input values arte PpasseGquas aera ecuce the ficw lcutemiee 
subreoutice CHKING, WhiGne1S Sed So 2a ore eee! a 
figures A.& throwgh (oo Uae eee ee 

After receiving «end verifying tne user's infut, the 
Frogram ktecins the calculaticns reéeyuired to deterfzire the 
cone sefi-vertex ancle. Frcozg tine o9 to dire 40C, tne 
Erogranr calculates the initial thermodynamic prozerties of 
tne water crI air, as tne case ay De, and calculates tke 
changes shich occur in tne fluid properties as a resi teem 
the fassage or the shock front. AS mentioned frevicusly, 
these calculations icr the air case are ratner straicnseee. 
ward and, kecause Of this Simplicity, tne computer CCaewaee 
the air calculation tas seen flaced iniine. For the »sater 
case, tke subroutines WSHOCK and EPV are needed to caictlate 
tne reguired thermodynamic properties. Sucrautine wWeooeee 
calculates the theraucdynamic fEreperties of the waver 
roint 1 upstream of the Shcck front and at foint 2, eee 
across tle shock frert on the downstream or shocked Side. 
Appendix A, figures 4.12 and A.13, contains tne £icweee 
which describes the Jlcegic of this subroutine. Subrcutine EPV 
was ccpied, with the euthor's permission, from Fuhs | Re toueem 
and is fully descricec in that work. However, toe FCR Dae 
traneélaticn of EPVY as included im the  srograms ii ae 
CON Welle ines e cmd xc. 

Faving calculated the initial thermodynamic tluid ¢GcnGae 
tions cn Ecth sides cf the sheck front, the grogran tCeaquie 
the iteraticn precess reguired to determine the cone seni- 
vertex argle. These calculations bpeyin at line 400 cf the 
computer ¢reogranm. This itératico process makes use cf the 
Eusegang graphical jategration technigue discussed Ey 
Sha pice eine em cae First, the radius of the Busecmarn cuewe 
for afpeint J is calculated uSing the “vetocr:, at “toa 


foint, tke streamline angle at that point, Samad tne wvaliewe- 
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Mitemorele Chedja dt that ~Olr:t {S€e Ligure 3.2). ext, fe 
oeormOl™@EreN—NClLCGlGerOl tac f£aSeaaik CUurVe at tre FoAant J 
mee G Ue atecqmictma tee velccity Comeonents at Ecint J 
feet ale Gesignacec U dna ¥}, the raiius Caicilated ir tae 


Erevicus step, ana tlre value cf the 


fo 

=) 
AL 

oa 


€ Omega. This circie 
M@iemmmer Lrcvides fhe fcimt FErCmM Which the pregran "dréwse" tne 
eeenwucsed tc calculate the mext Ecint On the Fusemann curve. 
Mmememyy tbe velocity; and the Susemann velocity ccaEcrenis, 
Mammemiexht POlnur, Jtl, are Calculated by "Swinying an arc" 
Meee the Ciltcle Center calculated in the previous ster. in 
Weaeertcn, the <stregiline angie at the next fEcint is 
@mereculated. 

The -~regram next tests tc see if tne cone surface ras 
keen reacked {line 475). This test is conducted Ly deter- 


Wininc if the arsoltte value cf the difference betweer the 


(b 
tae 


streazlire angle and the angle omega is 1ess thar a sfec 
i 


Ie 


| 


Mmrea test value (which is set to 1 X 10 78) . oe eh Sa 
ence is less tnan the test valve, the cone surface has ceen 
reached. In this case, the program then calculates tke ther- 
pwodynanic fluid conditions at the cone surface and tke ccne 
Semi-vertex angle ard displays the final results cif the 
EProgran. If tne difference Letween the streamline angle and 
the angle cmega 1S cleater than the test value, tne ccne 
Suriace Las not beer reached. Pieces c , tne Eregraa 
@erpoulates tHe thericdynamnic frivid properties at tae next 
Bernt, terete LOOpSs back tO  isine 400 to tesin tne 
meeeganh Graphical irtegraticn freocess for the new feint. Ia 
the case of a water run, tne subroutine WAIVEL is irvoxed 
Meee rricr tO leoping back tc line 400. 1p Si SUEECULL RS 
calculates the thermcdynamic frogerties of the water at any 
Meee; Girt. |Atrendis A, figures A.14 and A.~15, gives a 
Mmeemenalt whicn demcrstrates the logic ficw witnin this 


sme cutine. 
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ie 


Meme cr "“reloacging’ tke srograg Letween eaca executicn. 


\- 
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feature was included Ly solicitinjg a res,onse iroz the 


tn 
ry 


use 
at lire 15C€0 as tc whetner née/sne wisnes to make anctrer 


Mm 


L 
Peeenticn cr Ehe Srocram. An dadifirmative respcnse causes ta 
mmegrai t£¢ lcop crack tO Statezent 2t at tne Dbeginnirg cf tne 
Pmeelaiization Secticn of: the Main program. A negative 
Lespcersée causes [rogram termination. 

Nete that certair Separately compiled functions included 
mith this frogram have not Leen ftlowcharted because tneir 
Meme 15 £C StLlaigntiorward and because the sSudprocraagas in 
Mere tich usually consist of criy one or two lines of E€xecu- 
meee CCGG. The furctions fitting into this cCategcery are: 
nee, RICE, and SOSAIF. 


Pee ck LNSIRUCTIONS 


As disctssed in the objectives listed in chapter 1, cne 


Seas cf this thesis sas to ensure that the preogranx presente] 


was e€asy tc use and géintain. Ease of use has been facili- 
tateq ry including code which verifies that the inrut 
provided Ey the user 1S within specified ranges. eG 


Smmcie€, in the water Calculaticns, the pressure at rpcizt 

Gamers ce .creater tran 100.0 kilebars. This is due tc the 
Meme trac FrunsS' (Ret. 4] Sukroutines, which are usec to 
@mecilate tie thermcdynamic froperties “of the” water, are 
meeed Cn the wOrk OL Fat@hardsonyy et. al., [Ref. 3], which 
cnly gives results up to fressures of 100.0 kailckars. 
Tnerefcre, while tte subrcutines could compute the water 
SE@eaiticms at ~Ekessures greater tnan 100.0 xilobars, it is 
tieertain whether the values so calculated would ke entirely 
G@mrcct. FCrL this reaschn, the ranye of “the input values has 


teen restricted. Ease Om Wee is furtoer facilitated hy 


Ca) 
ta) 


CLOVEC ee weico e. Legh tgs ta eu & athe ASSGRCixX £ €CCMiceee 


(n 


fis 


sample cttprits generated oy tne ,.fogrtam@ £Cr the Suliags ee 
c 


comcléte ctcrint cpticrs Of “tne [eo lr aa eect GO eee eo 


tu 
es 
{2 


eee. 

Cne cf the awajor criticisss leveleaq sy computer Scene 
tists against the FCFIRAN language 1S Zor 1tS Lace 
reguirezent for formally defining all variarles usec witsin 
a preocragz. Similarly, FORTRAN haS an implicit declarauees 
policy satereby any variakbie whose first letter is I, J, kK, 


I, & cr N as imgrlicatly declared to be of Lyre Lrtege rae 


PD 


O 
rth 


cThLls Weeeearau, these two unfortunate characteristics 


}-~ 
fe? 


FORTRAN are avoided [Ey requiring explicit Gecladraticn One 
Frogran variables. 

Fase of Maintenarce and eéase oF understanding cf the 
Frogram tave been ackieved througn the use of a likerally 
commented frograz ard through the use of relatively feén- 
lngful variable names. in this regard, the FCRTEAN lJancuage 
is less thar desirable since it limits variacle names tc six 
characters. Appendix C contains a fully documented itera 
cr the ccmputer procram and all subroutines or LtLunctieae 
used Ey the main procram (cther than standard licrary furc- 
tions such as sir). Tables Vi, Vif, Vili and Ik, liccatcdiiee 
Appencix A, contain lists cf all variakle names used in Gme 
Taln f~rocram, thelr teaning, and their AkKS ({neter-kKilcocreme 
SeCouiais Gte.ts yer Byres see Sc Tacle X contains a Si 72a 
variarie list for the, Strroutine WSHOCK and Taric ie 
contains a variatie jist fer the Stbrobtine WNATVEDS ieee 
two takles ere also lJccated in Appendix A. 

Finally, this prcgram was developed using the FORTRAN IV 
danguace supplied by the LEM Ccrporation as part )ei tiem 
IBM 27CysQ232AP comptter syster which is the mMainmwech un 
systenr available at the Naval Fostyraduate School. As far as 
3S kncecwr, no impléegentaticn specific features have Leen 


inclugeq in eho Seema ame Thererore, tne prograr skculd 


34 


eeeteriee Goan Syste Woden Ha= a2 FSGazaati COLriler. icf c¢crcer 
Meme x ecco lte tae ELOGLah Om tne fo% 370 cemeut S 
eee eae tc C Nay a eet eatawe | SemoOl, @ tne TOLicwilts <sters 


Hust Le acccapliisneda ir the créer given: 
CDeene =o aeinah UUSem ko COtleried using the cotmanc: 
POrenAw@e CN EPL OW 


Meee CORHALG invekes the PORITRAN H Extended ccmeviler, wA1ich 
Mero rtih izing, “=~icducticon run Compiler supolied as part 
Gree tre IeM comruter system. The proyram could alse ce 
compiled using another FORTRAN compiler such as the FOFTEGI 
Semeediers. NOte that the abcve assumes the frogragz suprliec 
Pyethis thesis has keen entered into a file wnich has the 
filename CCNEFLOW anc whick is cr filetype rORTRAN. Nete 
further that this step need only be performed cre tine 


frovided errors do net occur. 


CZ) Next, the libraries of Standard subroutines and 
met icks supplied Dy the ccmputer center must Le attacted 
Pemmeecmt ite. Tmis tf accomelished by issuing the fcllcwirs 


Semmacc: 


GIOBAIL TXILIZ FCRIMOD2 “ACD2ZEH 


Memee tbat this step need cnly ce performed once [Erevided 


memertal terginaticns dc not cccur. 


ee oueenom: On PAN language, a c[cCrmatted inputycutyut 
Meo) statement reaqtires tne definition cf an IyoO device. 
Pemeeexaetp le, at gwost installaticns, oy default, I/C device 5 
meammusea TOL QM put “a@rd0 170 device 6 is used for cutrut. To 
ieeeeeze the Erogram ci this thesis properly, the fcllcwing 
comMarnds mtst be issued to define the I/0 devices tc Ee 


used: 


- 


FLILEGEiY Qos ties Cer fhe 
FILELEF O7 Tea¥ (225) 
As a result of the arcve CCitands, ail Cutlcit 2 iC cee 


1 


grroyraz will be sritten into a Zile on tre User" sat eee 
CONEFICk with a £iletype ci OUTPUT. ft is not  =tp ieee 
necessary te issue tle FILELTEr 96 colNteana shown abcvc) ae 
if net, all outgut sill be wiEittem @t the terminal anagem 
not ke availaple for brinting. ihe second file Gefi cies 


commard {FILEDEF 07) is abscluteiy recuired or the (~reguam 


will ©et Ceerate. This Lile definition tells the cci pues 
that JyG unit 07 is tte computer terminal. In the (kr @aemam 
cf this thesis, a@all input from the user 15 rejgucsted wie 


T/O writ Cy. Therefcre, if unit O7 1S not defined) iia 


ELOGLag cannot PeceivVegamy inpuE: 


(4) Fianally, the prcegrat can be executed by issuing 


the ccmuand: 
LCAE -COWN ES LOR Se 2 ew 


When procram executicn is completed, a cCiyptic nessacemee 
the fcrn Rk; T=C.017,6.01 16:43:00 will be displayed acu 
terminal. This is siazly a message from tae Cenputer eee 
ating systea indicating tnat the task just reguested nas 
teen ccuzpleted. The user can now utilize the system's Cl eme 
ating ccumands to review toe cutput Irom the executi coma 


the yma O fan 


if the user desires to execute the Frogram again, rettrn 
to steer (3) of the akLove procedure. If the gutisut Lica ae 
Erevicus executicno is ne longer needed, the same FILEDEF 06 
ccmmand can re issued. However, if the previous execution's 
cutgpet is reeded, the filetype of the FILEDEr 06 Ccnmeme 


should Fe changed (e.g. change OUTPUT in the above ccmmand 
to (GU cee ie 


eG 


IV. PRCGRAM CALCOLATION 3=SUITS 
ae EReGnar ene Slee 6£CR THE CALCULATION OF CCNiCaAl eel we 
AIF 


ome UecEcmee cela ee eemedleulates the LTiow Over a Ccche 
ope alt. Numerous prcgram executions were made icr the air 
case at various upstream Mach numbers ana for varicus snacck 
emere=5 ihe results chtained from these frogram executicas 
are summarized in Tatljies TI and ili. The numbers fredcuced cy 
tne fregrag were ccnyared for the variables snhowr ir the 
tables tc the taobles given an Kinney and Graham [fef. £], 
Big@eamenet the yrarphs civen Dy “@hapiromj Ref. 8], ard tc the 
memes Eroduced by keral [f ket. 10]. Ine values given in 
these scurces are  a4ircluded in Tabies II and III, aceEre 
meerectriate, aleng with the calcuiated results. As cén be 
seen fren toese two tables, Mmmemedlemlations tade Cy tbe 
Frograr cive yuite accurate results (less than a 1% €rrer in 
Meee cases).  Certairly, the procram is much more accurate 
mmo cHoe'S ability tc pead the graghs presented in Shapiro 
feet. Gt. 

PoaoeOnGGmemeceoncChztakisens, 126 is believed tnoat the part 


JS 


Gene P~lCOman Which calculates Supersonic itlow over a ccne 
Mew ne=  aGewrate. Ihe importance OF this ract is tnat, Ly 
Peeling contident that the procedure followed for tke air 
case is valid and that this procedure has been ccrrectly 
jaoplerented in the froygrarming language, it 1s safe to 
assume tkat the Eusegann calculaticn procedure utilized fcr 
tne air case can ke accurately applied to the water case 
Meovided the water  Ccnditicns at each point are calculated 
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E. EFOGEAS nESULTS FCR GHe CS2lCuLragion CP ecg eee 
RATER 


As @fentionped 12 ne ~Pevileto Ss -eriew the Calculcavige 
Frocecure was verifiec for the air case EY CoOnmecalri=ca ee 
known results. This ccmMparison showed the metnod usec ir the 
calculaticns was correct. The next ster taken was tc deter- 
Tine if the subroutines used to caiculate the tnermcayranic 
Fropertiés cf the water gave results CChpararte tc tinccouee 
Racharasei ect. . Jae [Beta and Fuhs [ket . ge 
Calculaticns were jperformec which gave resultS which could 
Eke ccupared to Table Ii of Funes [Ref. 4]. Table IVY Eresenae 
a ccmyariscr betaeen the numeric results calculated fry tke 
Frograg and those given by Fuhs. aAsS can be seen, these 
results are, within reasonakle error, remarkably similat. 
This is nct unexpected since the subroutines used fy the 
Tain froc¢raz to calculate these varlarles dare ,esser riage 
direct language translations of Fuhs’ programs. Teererauee 
the results should be similar. 

Havirg verified thnat the procedure utilized in «the 
calcujaticn of the cone sSeni-vertex angle was ccrrect 
(thrcugh the aixr™ results) and having veriried that tae 
rrojram correctly calculated the therwodynaric J;rcrertie= ee 
the water at any pcint, it Was Deliliewed that the ~receae 
could ke executed fcr the warer case, FOL’ var wous Wa eee 
conditicre, with cértdainty that the results So caveuwi age 
would re accurate. Ecwever, a question arose as tc whetker 
the calctlations perfcrmed in the water case were indeéepen- 
dent cf tne upstream thermodynaric properties of the water. 
In adage tke only gtantities needed to calculate the ccne 
Semi-vertex angle are the upstream Mach numker (4) and the 
shock angle (oo ) (€.J- see figure 17.744) of ~ Shacagee 
peers 63am Linus, 6 Lot scum the Calemiaticn ef tie Cem 
semi-vertex angie 1S independent of the upstream temperature 


Cr pressure, and, therefore, universal curves can be drawn 


4C 
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which “describe tne cone Flex Under eee... ee 
enown Bye sSktapiro [ Réficas |; anetner universal results for 
the water case could ce O£tainNe€O was Rot sco cicac.: a tS 


ECiDLE if few Gisen ssc ee 

Li€grarn and rostko jRef. 11) discuss the general tner- 
nodynarics of fluids and introduce the concert Gram 
"“thergal" ecuaticn cf state. In jeneral, the therial eotee 


tion cf estate fcr any fluldei cee = 
E(p,0,1) = 9 (4.1) 


For a@ perfect gas, the thermal equaticn of state is given 


Cs 

DY = en 4 (ay, Ze 
where & TAS Anotréer equaticn which relates the Vtherieee 
Lamic variarles é, Vv, and Tis the "calorie" eu vatienmee 
state, which, in a géneral form, 15 given Ly: 

£ (egy, 1), eae (4.3) 


A Caicrically perfect gds 1s cetinecas. 2 


e = cil (4.4) 


f= a ee 
SS ( } 


where € if internal erergy and h is enthalrfy. Thus, the neat 
capacity (either Cc. or Sad 1s 4° Gomstame Rome cake: 1 came 


Ferfect cas. 


-~ 


P2e,cgamm ana temo ~{Rer. 11] wUuemenstrate trat tne 
[moreeies é¢L¢d Ca loOrll Gamecua a zm 


Mmerated Ly: 


(4.6) 


Mmemeeccucing equaticno 4.2 itte equaticn 4.6 results ir the 
Mmet=oarad Side of equaticn 4.6 becoming ¢€yual tc <2zeéro. 
Wms per crl a thertally crerfect cas, tne enthalpy 1S iadepen- 
dent cf the pressure, and, femceumecuGialoy 15 a Lurcticn 
Oeeey Gi tenrerature. Therefcre, Ompeecesscak Jecoma 1 ticE Lor 
Sumatacnh 4.5 to bre valid is to aave a thermally perfect Jas. 

memewaren, Richardson, et. dai.,,[ Ref. 3],. demcnstrate 
that tke Léeat cafacities are functions of the temperature cr 
the water. Thus, ecuations 4.4 and 4.5 are not valid fcr 
Meek, and walter 1S fet a Caloricaliy perfect fiuid. 

Ifa fiuid is Ecth thermally and caiorically ferfect, 
MmmemecnaiticnS acrcs< a norgczal sheck rrent depend orly or x 
= C)/¢, and on the freestream Mach nuaker. Conversely, if 
Miemtetic 31S not thermally ard Caloricaily perfect, addi- 
Ionia Vdaliablies must be specified in order to déeiire the 
mum@iax CCLGItTZOnS. Extending the argument, cre can state trat 
Meet pcLochac CLlOw Ck a therfaliy and calorically ferrect 
ieeea OFEr a Cone if a fuUrcticn only cf the heat capacity 
mec, Kk: the greestrean Mach number, 4G andar t 6 €. sSee « 
feeomte ancle, Gg. AS a result, a Single graph 15 sufficiert to 
represent all flews cver the cone. This is the conditicr icr 
empoecriect gas. 

Papeete wie mmeace rr LUlde1s Mat thermaliy ara calori- 
cally perfect, tne gclution £cr Supersonic conical flew is 


gependent upon variables cther taan Mk, M and.2 . 


l f 
Consecuenrtly, auniversal graph of the shock angle versus 


the freestream Mach rcurker, hith the cone angle as an entry 


f= 
(a) 


rarazreter, 20€S not €x1St. Sanee Watem 15 Hevea wee cee 
and calcricall; geriect fluid, it 1S 9n0t Qeossilt ic toma 
universéel ¢crapes OL UlLe, Van tas es eo eee 

I¢ gain insight into Eke SENSItTLIVIity of the SUEeCE See 
conical flew in watezr to the upstream conditions, 4 saa 
was ccndtcted in whick caiculatlons were tericimed tc deme 
Eine soetker the cone semi-vertex angie is relatively inde- 
Fendert cf the upstream temperature or whether vatricus 
Erogram executions are required tO SnOwW the Varilati comes 
cone sé€ni-vertex angle with water temperature. 

Tre cétermtinatior discussed above was made Eby helding 
the upstrear pressure, the shock angle, and ts¢€ urstreag 
Fach “DUIser COoOneecam e- Then the upstream temperature .as 
varlec tc céternzine the variation in tne ccne semi-vertex 
angle. In crder to Ecld M, ccastant, the Pressure [uae 
caticn factcr, described in Chapter I1.D, was variea Simao 
the Mach Runmber upstream for the new temperature fatcned 
(withir reasonacle accufacy) the upstream Mach tual e reas 
the crigina] temperature. Note tnat "reasonable accuracy" 
for gatching of the upstreaa Mach numbers meant Matching the 
hugcers to the fourth decinal fFlace. Table V prevides the 
PesuItS ire kh tivese acatecu lace co=. AS can be seen [r¢i che 
table, the core sé€mi-vertex anyles calculated fcr _ the 
Varicus tenperatures are all cf approximately the Sarre value 
(withir 124jy. The fressure dcwastream of the SaCCK eee 
varies Fy 153% as can fe seen in Table V. While there is small 
variaticr keéetween tke values for the cone semi-vertex 
anyles, it must Ee renembered that, for practical j;urpeseee 
these variations are very szalljl and can be safely ignecred ina 
€xamlininc the flow cver a cone in water. Thus, 1t was 
concludeda tkat the calculaticn cf the cone semi-vertex angle 
in water can be conducted inderendently of the water teuper- 
ature | (armieeeis it eee i the accuracy reguired for the cai culae 


tion cf tne moverent cf the metal jet from a shaped-charge). 
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thnrouch tre mecnhaniei cl ENE EieSseute Bll ewe scat epee 
Eavinc determined that ta€ Caicuiatiorns Cat Ss] cereus 
indejyendenctly of the temperature, plots of the variakles of 
interest in water s€re made in a manner similar te the 
cracns fc~L @1©& ShOWD IU @shatelre kee oe - 

Finaliy, hnaving cetermined that the calculaticns couse 
ke mace independently ci the upstream temperature ard t£at 
these caiculations sculd be accurate {within 1% fer cche 
angle and 12% for pressure), i1t was decided to pict varicus 
Farazeters cf irterest for the flow over a cone in sater. 


These créphs are presented as fcllows: 


(1) Figure 4.1 is a plot cf the £reestream Mach Nungter 
versus tne shock angie with the cone semi-vertex 
angle as an entry rarameter. LThis  grath Ls eee 


water analog to Figure 17.7.{a) of Shapiro [| ket Seage 


(z) Figure 4.2 is a plot cf the freestream Mach rugter 
versus the Mach numker at the cone surface with the 
cone semi-vertex angie aS an entry parameter. Tris 
Clagph 1S (‘the water™ analog uo Seeger ei nee OL 


Sieyetes Pees 64). 


(2:) Finéelly, Figrre 4.3 is a Elot of Ehe Ereest recneiaes 
DUMEEK VELSUS Che Ciao Cecio meres ee witkh the 
ccnhoée semi-vertex angle aS an entry parameter. This 
¢craph 1s the water analog to Figure (7257 = (ae of 
Shapfiso [Resé . sane 


It shcujd be nected here that tne Graphs deveicred “ama 
Fresented in this tktesis were obtained through rerpeated 
executicns cf the CONEFLOW fFrogranm, which Generated sa pros. 
Rately 3€CG datd Seintismecr semen. these data (pot mre. wei 


Fre-frccessed -te rencve entries with a Mach nuster oF 


4E 


Meeucimesen  O.0 {2 cCenveriegyce £O> case CL  Presertéatior 
ire) . ie Se eeG@ec sO oeemmec ta) | Ol LtS Yores cassec tc ar 
Pete eC lctrecauy aha. STOCEnINS SEOgGran csvyelo-ea cy fraike 
Meee |Z nner ised a@ GNadratic Snhepara’s metnod tc sascta 


tt 
O) 


‘ 


Oy 
ih 


mpe cata ECintS generated’ After soothing, tne da 
post-[iIccessed to remove erroneous iainterpelaticnr fecirts 
wnich were senerated btecause the original data sas irreg- 
Glar. Fanaliy, the ~rocessed data was passed to a FCFTSAN 
prograz developed Fy the author or this thesis which 
mmmieaed §~the COLtcUuring feattre or the DISSPLA grarring 
Meme ares Or the ISSCC company. ime results of tnis extensive 
Erocessirg are the graphs presented in this thesis. eE€catse 


Memeetse Itregularity of the criginal data {which causes 


@imciculties in tre Meetid cima Lac Suiepard's intersciatica 
methcd), and because of the smoothing reguired Ey the 
@emeouring routines cf DISSeLAa, Bie wGbatasS presentec are 


Somewhat rcugo and skculd, tneretore, only be used rer "cacx 
c£ tke ervelope" calculaticns. Accurate Cresuits are 
mmeyvrced Ly the CONEFLOW pregram which snould be utilizei 
meme CLE precise work. 

mindiily, Ole OL the primary reasons Lor translating the 
Mmeeeitel EPeasic Ctrogrems into @ higher-order languase {in 


additicr te the desire to make the programs nore acceéssi} 


ty 
pas 
48) 


Meroe Steed the execution ci the progran. iaeceses re ca 
Meme work Cr this tresis has more than accomplished tnis 
Merete ft MUSt bE Menticned sere that tae ELregranr will 
execute extremely iast, considering all the iteraticns 
needed, previded the computer system in use is net heavily 
in use at tke tine. For example, at the Naval Postgraduate 
Sschocl, when the ccnrputer system is relatively free, the 
Frogram cperates so ctickly that one execution requires less 
meas seconds {irom the time cr the last user inpgut te the 
Meme €t the Fequest by the fErcgram for an indicaticn of 


whether ancther execttion is desired). Conversely, when the 
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Wa) CGN ereUreene nies 

As discusSseq if Schein rcdue tren tnis thesis aad 
certair scalsS whick it was desired to achieve. it ime 
kelieved that these coals have ail been successrully met. 
The pregrar of this thesis has EFeen successfully trarslated 
from Easic into FORTEAN as desired. In the translatica of 
the criginal preygfan, modern structured pregranmging eae 
tices have been folicwed tc the greatest extent pessitle. 
Finally, tke ~rogragy presented 1s guite "wser-rriend] yaiags 
is very well-documented. 

As was demonstrated in the resultsS chapter, this [Ercgqiag 
calculates’ correct results icr €ach cf the air Cr ieee 
cases. It iS mertioned here that, as pointed out frequentl; 
Ey ccemputer scientists, it is virtuaily ifposSsibie tc i quam 
rrograg fcr all possikle cases. Therefore, no f[rCgranwe 
entirely error-free. It is kelie€ved that the prograt ci trae 
thesis is as free of Errors as 1S possible without exrnaus- 
tive, tize-consubing and extremely expensive testing. 

As a resuit of ELaving wet tne three goais menticned in 
-the first fLaragraph, it is Leiieved that the progtamt Gite 
thesis Wali) prevLac¢ an excellent WOrTKAING tcGu ioe 
researckers in the field of shaped-charye jet penetraticzc int 
water. The fact tnat the pregram is well-documented wae 
take any medificaticrs or extensicns to the progranz,) srouie 
that ke required or desired, easy to perform. Further, sites 
the pregrat is €asy to usé, cnly acursory kuacwledge of 
computers is needed in order to utilize the frogran. tMIhese 
featules aLe TequiregentSs LCr any COMpUeet programe.) 1 core 
to be used wey Sciemtists astarytool. 100 or ten ece tetera 
their pregrams require that the people whe simply desire to 


use their capabilities must lé€arm a Great deal “ci detaw 
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meemdive dLgeundernstanding cf FORIRAN 


Mae oe Zech Lnhe  cbOGno me sC@mtes Gesigcn. Uritcrtunately, 


Ww 


Miiemw tc tle Crebaring system cr the i314 3/0 Computer Svster, 
mies sdfie Cannot re Said about the stess reguired in crdeér to 
fememaily use the Gpregrat. Scme familiarity witn the ccrguter 
system in use at the users Iccation will cre a necessity in 
mmmer tc cpérate the f~regragi ccrrecti;. 

Pali, Ao echt LoOmegu ln Lae “—“intreduction tc tris 
thesis, tke program sili previde an excellent test case ani 
@emeatisch wodel for a COMEUTEC program vwhicn mcdels the 
peered cLicw OL tne metal jet irom an explosive shaped-ckarge 
fired thrctch the water. Picmacmual Siltuaticn 1S a Elurt= 
Mosed, Latker than 3 Sshare-~cinted conical, ier ee eee rn 
mipben 12 Euco more difficult tc solve. Therefore, a knecwr 
Bemuticnh and méethodcicgy or sclution wfor the easier rrcklen 
mea Necessary first step to the solution of the larcer 
Meeker. It 1s believed that the program of this thesis will 
serve as this necessary first step. Further, it is relieved 
epat, wken the equaticns solving the actual flow frorkleg are 
develcred in their final fora, ENememewEan Of this thesas, 
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